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In a study of B + — » J/tp^K + decays, we find evidence for the radiative decay A"(3872) — » J/07 
with a statistical significance of 3.4cr. We measure the product of branching fractions B(B + — > 
X(3872)^+)-B(X(3872) -> J/07) = (3.3 ± 1.0±0.3) x 10~ 6 , where the uncertainties are statistical 
and systematic, respectively. We also measure the branching fraction B(B + — » = (4.9 ± 

0.2 ± 0.4) x 10 -4 . These results are obtained from (287 ± 3) million BB decays collected at the 
T(4S) resonance with the BABAR detector at the PEP-II B Factory at SLAC. 

PACS numbers: 14.40.Gx, 13.20.Gd, 13.25.Gv 



The X(3872) state was discovered by the Belle Collab- 
oration in the decay B+ -> ^(3872)^+ 2]. This 
signal was confirmed by the BABAR Collaboration 
as well as the CDF and D0 Collaborations Q. Inter- 
preting this new state has been challenging. Its narrow 
width, mass near the D°D*° threshold, and small branch- 
ing fraction for the radiative decay X(3872) — > Xci7 have 
made it difficult to identify the AT (3872) with any of the 
predicted charmonium states 0|. Alternate proposals 
have been made, including a D°D*° molecule 0, or a 
diquark-antidiquark state |7J . Evidence for the radiative 
X(3872) -> decay in B+ -> X(3872)K+ would 

determine the C-parity of the X(3872) state to be pos- 
itive, limiting the conventional charmonium assignment 
options while remaining consistent with D°D*° molecule 
model predictions. 

A number of other new states have recently been 
found. The Belle Collaboration has claimed the discov- 
ery of a broad resonance in B decays, referred to here 
as the F(3940) state 8|. The nature of this state is un- 
known, and there is no reason to yet preclude B + — > 
F(3940)iC + , F(3940) -> J/07 as a possible decay chan- 
nel. Belle has also identified a possible x' C 2 charmonium 
candidate in two photon production, referred to here as 
the Z(3930) state 0. This state could be produced in B 
decays, and if the tentative x' C 2 charmonium assignment 
holds true, it should decay radiatively to J/07 (albeit at 
a rate predicted to be quite small). 

We study the decay chain B + — > ccK + , where cc de- 
cays radiatively to J/07, and the J/0 subsequently de- 
cays to a lepton pair. The notation cc represents conven- 
tional charmonium, such as the triplet % c j(lP) states, 
or any state with positive C-parity for which the J/tpj 
decay is not forbidden. 

The data sample for this analysis consists of (287 ± 3) 
million BB pairs collected with the BABAR detector at 
the PEP-II asymmetric e + e~ collider. This represents 
260 fb -1 of data taken at the Y(AS) resonance. The 
BABAR detector is described in detail elsewhere ^] . The 
innermost component of the detector is a double-sided 
five-layer silicon vertex tracker (SVT) for precise recon- 
struction of B-decay vertices. A 40-layer drift chamber 
(DCH) measures charged-particle momentum. A ring- 
imaging detector of internally reflected Cherenkov ra- 
diation (DIRC) is used for particle identification. En- 
ergy deposited by electrons and photons is measured 
by a CsI(Tl) crystal electromagnetic calorimeter (EMC). 



These detector subsystems are surrounded by a solenoid 
producing a 1.5-T magnetic field. The flux return for 
the magnet is instrumented with a muon detection sys- 
tem composed of resistive plate chambers (RPC). For the 
most recent 51 fb -1 of data, a portion of the muon sys- 
tem has been replaced by limited streamer tubes (LST) 

A J/0 — > £ + £~~ candidate is reconstructed by combin- 
ing a pair of oppositely charged muon or electron can- 
didates having an invariant mass compatible with the 
nominal J/0 mass. An attempt is made to recover en- 
ergy loss due to bremsstrahlung by searching for photons 
near electron candidates. Candidates for J/0 are then 
combined with a candidate kaon and a photon to form a 
B+ — * J/ip-fK + candidate. 

The J/0 — > e + e~ candidates are formed with elec- 
trons (and bremsstrahlung photons) with 2.950 < 
m(e+e"(7)) < 3.170 GeV/c 2 . Candidates for J/0 -> 
require muons with 3.060 < m(/i + /i~) < 3.135 
GeV/c 2 . The cc candidate is reconstructed from the 
mass-constrained J/0 and a photon with Ey greater than 
30 MeV. Additional selection criteria are applied to 
the shape of the lateral distribution (LAT) |l4J and az- 
imuthal asymmetry (measured by the Zernike moment, 
A42) 0| of the photon-shower energy deposited in the 
EMC. The radiative 7 candidate is rejected if, when com- 
bined with any other 7 from the event, the invariant mass 
is consistent with the 7r° mass (see Table Q. The ratio of 
the second and zeroth Fox- Wolfram moments (R2) 0] 
is used to separate isotropic B + events from typically 
anisotropic continuum background events. The mass of 
the cc candidates, m C c, is calculated by constraining the 
B + candidate to the nominal B + mass. 

To identify B candidates, we use two kinematic vari- 
ables, m,B and m m i ss . The unconstrained mass of the 
reconstructed B candidate raj = yj E\ /c A — p 2 B , where 
Eb and ps are obtained by summing the energies and 
momentum of the particles in the candidate B meson, re- 
spectively. The missing mass is defined through m m i SS = 
\/ {Pe+e- ~ Pb) 2 , where p e + e - is the four-momentum of 
the beam e + e _ system and Pb is the four-momentum 
of the B candidate after applying a B + mass constraint. 
These variables are uncorrelated by construction, and are 
advantageous for analyzing B decays in which a particle 
in the final state has poorly measured energy. Events 
with a correctly reconstructed B + decay should have val- 
ues equal to the nominal B + mass for both kinematic 
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TABLE I: Summary of acceptance criteria for candidate 
events. 

Region 

2.950 < m e+e - (7) < 3.170 GeV/c 2 
3.060 < m„ +M - < 3.135 GeV/c 2 
R 2 < 0.35 
LAT < 0.5 
A 42 < 0.1 
Reject 122 < m n o_^ < 145 MeV/c 2 



variables. 

To best separate signal from background, the signal 
selection criteria are chosen based on Monte Carlo (MC) 
samples by maximizing the figure of merit ns/(a/2 + 
y'rag) [l7| where ns and ns are numbers of signal and 
background events, respectively, and a represents the 
minimum level of significance desired. For this analy- 
sis, a = 3 is chosen. The optimization is performed by 
varying the selection values for m e + e -( 7 ), m^+j,-, i?2, 
photon LAT , photon A42, and the photon 7r° veto, while 
requiring tub and m m i SS to be within 100MeV/c 2 of the 
nominal B + meson mass. The optimized criteria used in 
this analysis are summarized in Table [I] 

We extract the signal with an unbinned two- 
dimensional extended maximum-likelihood (ML) fit to 
the kinematic variables tub and m m i SS in 10 MeV/c 2 bins 
of m C c. Fits failing to converge or lacking statistics 
are combined with adjacent m C c bins to ensure fit suc- 
cess. The probability density functions (PDFs) for sig- 
nal extraction are the product of independent fits in tub 
and rrimiss , defined separately for signal and background 
events. 

The signal PDFs are determined from Monte Carlo 
simulation of B+ -> Xc\K+ and B+ -> X(3872)JC+ de- 
cays. The to b and m m i SS distributions of B + — > ccK + 
signal events are both modeled by a functional form 
similar to a Gaussian with asymmetric tails, /(x) = 
exp[— (x — m) 2 / (2u 2 _ + a±(x — to) 2 )] , where the ± sub- 
script indicates different parameter values on cither side 
of the central peak. The signal PDFs for these two cc 
modes are found to be equivalent to one another within 
statistical uncertainty. 

The background consists of two parts, a combina- 
toric component with a flat distribution in the kinematic 
variables ms and rrimiss, and a component that peaks 
in m m iss composed of B decays similar to our decay 
mode. The peaking background events are mostly from 
B+ -> J/ipK+ir , tt° -> 77 and B+ -> J/tpK*+, K*+ -> 
if + 7r°,7r° — > 77 decays. These events are incorrectly 
reconstructed as the desired final state if one of the pho- 
tons from the tt° decay is undetected. This background 
does not peak in the other kinematic variable ms, nor 
in m C c- The only doubly-peaking background may arise 
from B + — > JipK* + ,K* + — > K + j. These events can 



peak in both ms and m m iss, but the branching fraction 
for this decay mode is small and can be ignored. The sim- 
ulation also indicates that the combinatorial background 
is almost entirely due to B decays. 

The background PDFs are fitted to events from generic 
B + B~, B°B°, qq (q = u,d,s,c), and t+t~ MC samples. 
In rriB, all background events are modeled by the tail 
of a wide Gaussian function. The m m iss distribution of 
background events is parametrized by the ARGUS back- 
ground shape (lsf for the combinatoric component, while 
the peaking component is characterized by a Gaussian 
function. 

The maximum-likelihood fit returns the number of 
B + — ► ccK + signal events, N S i g , in each 10 MeV/c 2 m C c 
bin. The number of signal events is found by fitting to 
the N S i g versus m C c results with functions representing 
the cc mass distribution of each signal mode. Based on 
Monte Carlo simulation of the Xcj=o,i,2 

and AT (3872) [jjj 

decays, the to cC shape for each of these signals is individ- 
ually parameterized with a double Gaussian distribution. 
In the fit to the ML results, the Gaussian means, widths, 
and ratios of the areas are fixed to the values determined 
from the MC simulation, with the heights of the peaks 
permitted to float. As N S i g can also include non-signal 
events peaking in both to b and m m i SS , a first order poly- 
nomial in N s ig versus m C c was included to account for 
the level doubly-peaking backgrounds. The number of cc 
events is calculated from the area of the fitted Gaussians 
above this background. 

The effectiveness of the signal extraction method is val- 
idated on Monte Carlo samples for Xco.i and X(3872). It 
is found that the proximity of the large Xci signal peak 
introduces a significant bias in the measurement of a Xci 
signal with this method. Therefore we do not quote re- 
sults for the Xci mode. Successful performance of the 
AT (3872) extraction is verified on Monte Carlo generated 
samples for numbers of events similar to the measured 
value, as well as for the case of a null result. 

The efficiency is determined by calculating the fraction 
of the events generated in Monte Carlo simulation that 
survived the analysis selection criteria from Tableland 
are returned by the fitting procedure. Standard BABAR 
corrections are applied to account for particle identifica- 
tion and tracking differences found between simulation 
and data. These corrections are at the level of a few per- 
cent. The resulting efficiencies are (16.8 ± 0.2)% for the 
AT(3872) mode, (13.3 ± 0.2)% for x c o and (13.5 ± 0.3)% 
for Xci , where the errors are statistical. 

Systematic uncertainties on the branching fractions 
are reported in Table [H] Sources include uncertainty 
in the number of BB pairs, uncertainty in the secondary 
branching fractions for Xco.i Jfyl and J/if) — > £ + £~, 
PDF parameterization uncertainty due to MC statistics, 
uncertainty in the m C c parameterization, particle identifi- 
cation, tracking and photon corrections, effects due to fit 
technique (such as choice of m C c bin width and fit start- 
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TABLE II: Summary of systematic uncertainties. The uncer- 
tainty due to secondary branching fractions (BFs) does not 
apply to the product of branching fraction results. 



Source 


Xco(%) X 


ci(%) A(3872)(%) 


B counting 


1.1 


1.1 


1.1 


Secondary BFs 


8.5 


5.4 


1.0 


MC statistics 


16.5 


3.2 


8.7 


m C c shape 


3.1 


1.3 


1.5 


Particle ID 


2.0 


2.0 


2.0 


Tracking 


3.6 


3.6 


3.6 


Photons 


1.8 


1.8 


1.8 


Fit technique 


1.7 


1.7 


1.7 


A(3872) mass/width 






2.0 


Total 


19.5 


8.1 


10.3 



-10- 




3.7 



3.8 



3.9 



4 4.1 
mfGeV/c 2 ) 




(GeV/cr) 



FIG. 1: Number of extracted signal events versus m C c for Xcj- 
The solid curve is the fit to the data. The \ 2 P er degree of 
freedom for this fit is 48.7/34. 



ing point), and uncertainty in the true X(3872) mass and 
width. The uncertainties due to MC statistics, m C c and 
AT (3872) mass were evaluated by varying the individual 
parameter values by la from their measured values and 
finding their effect on the signal yield. The largest source 
of uncertainty (aside from secondary branching fractions 
which are beyond the control of this analysis) is due to 
the variation in signal yield with the choice of PDF pa- 
rameter values. In the case of the X(3872) signal, the 
total uncertainty is dominated by statistical rather than 
systematic errors. 

Figure ^ shows the fit to m C c in the mass range 
3.311 < m c s < 3.711 GeV/c 2 . We extract 27.9 ± 11.7 
XcO events and 807.2 ± 33.3 Xci events. Using our sig- 
nal extraction efficiencies, we calculate the product of 
branching fractions B(B + — > Xc\K + ) ■ B(xd — ► J/07) — 
(1.76±0.07±0.12) x 1CT 4 and B{B+ -> XcoK+) -B{x c a -> 
J/07) = (6.1 ± 2.6 ± 1.1) x 1CP 6 , where the first error 
is statistical and the second is systematic. Taking the 
branching fractions for Xc0,i — > J/07 from [l3|. we cal- 



FIG. 2: Number of extracted signal events versus m c - for the 
X(3872) region. The solid curve is the fit to the data. The 
X 2 per degree of freedom for this fit is 57.8/37. 



culate B(B+ -> xaK+) = (4.9 ± 0.2 ± 0.4) x 10~ 4 , 
and B(B+ -> XcoK+) = (4.7 ± 2.0 ± 0.9) x 10~ 4 cor- 
responding to the 90% confidence level upper limit of 
B(B+ -► XcoK+) < 7.5 x 10~ 4 . These branching frac- 
tion results are consistent with the current world average 
[l3j . and in the case of B + — > Xc\K + , more precise. 

We extract the number of A" (3872) signal events in the 
mass range 3.672 < m C c < 4.072 GeV/c 2 and find 19.2 ± 
5.7 events (Fig. We derive the product of branching 
fractions B{B -> X(3872)K+) ■ B{X{3872) -> J/07) = 
(3.3 ± 1.0 ± 0.3) x 10~ 6 . The statistical significance of 
this signal, taken to be the square root of the difference 
in x 2 values between the fit in Fig. |3 and a similar fit 
assuming zero signal events, is 3.4er. 

Additional fits are performed to search for the F(3940) 
and Z(3930) states by adding a resonance in the appro- 
priate mass region. The measurement of the F(3940) 
state from Q finds a mass of 3943 ± 17 MeV/ c 2 and width 
of 87 ± 34 MeV/c 2 , while the Z(3930) state is found to 
have a mass of 3929 ± 5 MeV/c 2 and width of 29 ± 10 
MeV/c 2 0, where the statistical and systematic uncer- 
tainties have been combined in quadrature. We model 
the mass resolution for the decays of each of these states 
to J/07 by a Gaussian function in m C c with the mean 
and sigma fixed to the Belle measurements. Because 
the masses and photon energies are similar, we assume 
the same efficiency for these modes as for the AT (3872) 
state. We find — 16 ± 34 events and —5.4 ± 8.3 events 
for the F(3940) and Z(3930) states, respectively. We 
define an upper limit on the product of branching frac- 
tions by assuming a Gaussian distribution for the number 
of signal events and its uncertainty, and integrate over 
the physically-allowed region from to 90% of the total 
area around the mean. Systematic errors are estimated 
from the contributions listed for the A"(3872) in Table 1TT1 
Uncertainties on the F(3940) and Z(3930) masses and 
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widths dominate entirely. The total systematic uncer- 
tainty on the product of branching fractions is 101% for 
the K(3940) and 22% for the Z(3930). To account for 
the width uncertainty, it was varied by la from the mea- 
sured value and the largest resulting upper limit retained. 
Using these basic assumptions, we calculate B(B — > 
y(3940)X+) • £(T(3940) -> J/0 7) < 1.4 x 1CT 5 and 
B(B -> Z(3930)if+) • Z?(Z(3930) -> J/^7) < 2.5 x 10~ 6 
at the 90% confidence level. 

In summary, we measure the branching fraction 
B(B+ -> Xci^ + ) = (4.9 ± 0.2 ± 0.4) x 10" 4 and de- 
termine a 90% confidence level upper limit of B(B + — > 
XcqK + ) < 7.5 x 10~ 4 . We find the product of branching 
fractions B(B -> X(3872)K+) ■ B(X(3872) -> J/tpj) = 
(3.3 ± 1.0 ± 0.3) x 10~ 6 , with a statistical significance 
of 3.4a. This provides evidence of the radiative decay 
A (3872) — > J/^7 and of charge parity C = + for the 
A (3872) state. We search for radiative decays of the 
Y(3940) and Z(3930) states to J/ip 7 in the B+ -> ccK^ 
channel, and find no evidence for such modes. 
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